Inhomogeneous states and nodal fermions in the SU{2) gauge theory 
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We discuss the issue of phase separation in the SU (2) slave boson theory of Wen and Lee of the doped Mott insulatop^. It is 

shown that the constraint structure of the theory leads to the interpretation of the holons to have hard-core interactions, 
which is demonstrated further by studying the empty limit (no electrons). Surprisingly, with hard-core interactions even the 

empty limit is described well by the slave-boson theory, both as an energy density and with the regard to dynamical 
properties. The consequences are investigated in the overdoped superconducting regime, where both phase separation and a 
d-\- s structure of the order parameter is obtained. This s-wave component is already imminent in the description of the hole 
in the slave boson theory. The interacting nature of the holons also lead to sound modes in the single-electron propagator. 
The novel idea of the isospin spiral is introduced, based on the projective symmetry principles of Wen. This isospin spiral 
explains the coexistence of superconductivity and the Mott insulating state, being the consequence of phase separation. 
Secondly, it might be able to explain why nodal fermions survive in the presence of charge inhomogeneities. 
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I. INTRODUCTION 

A long standing question in the physics of high- Tc su- 
perconductivity is how nodal fermions can coexist with 
stripes, as experiments seem to point out^. There is no 
simple theoretical explanation on the market as to why 
charge inhomogeneities do not interact strongly with the 
low-lying degrees of freedom. A possible explanation is 
that in the high-Tc's, electrons are splintered into spinous 
and holons. Being different degrees of freedom, it is pos- 
sible that charge inhomogeneities do not communicate 
with the fermionic low-lying excitations. Experimental 
support for such fractionalisation is for example given by 
photo-emission- . 

This idea of fractionalisation is incorporated in the 
slave boson theories, leading to spin liquid states^i^. 
These spin liquids are states "nearby" the parent Mott 
insulating states, and carry nodal excitations'^. The hope 
is that these spin liquid states are candidates to describe 
the evolution from the Mott insulator to the d-wave su- 
perconductor, in which quasiparticles with d-wave dis- 
persion are extensively proven experimentally^. 

These slave theories are plagued by one problem, how- 
ever: they appear not to be able to describe inhomoge- 
neous states, like stripes, let alone that they can explain 
coexistence of nodal fermions with stripes. The main 
point of this article is that the SU{2) slave boson theory 
as developed by Lee and Wen^i^ii^, is able to capture both 
aspects. This conclusion is based both on the projective 
symmetry ideas of X.-G. Wen^°, and on an improvement 
of the original theory by Wen and Lee. 

The first idea is at the root of our proposal of 
the isospin spiral. Wen showed that for zero doping, 
the mean field states describing a staggered flux phase 
(SFP)^^ and the c?-wave superconductor (dSC), are phys- 
ically the same, i.e., they are gauge equivalent. Only for 
non-zero doping, this equivalence is broken, such that a 
d-wave state is favoured, mimicking the instability of the 
SFP against the dSGi^. Our idea is now that for low 
doping, an inhomogeneous state is favourable above a 



homogeneous dSC. This state connects an SFP with no 
doping smoothly with a charged superconductor. This 
leads to a picture in which isolating SFP states coexist 
with dSC states on stripes. The protection of the nodal 
fermions lies in the fact that for zero doping, both the 
SFP and dSC carries nodal excitations. 

For these inhomogeneous states to exist, it needs to 
be proven that the SU{2) slave boson theory supports 
phase separation, which is not proven so far. We show 
that this phase separation does occur, which is con- 
nected to our technical improvement of the original Wen- 
Lee theory. It is demonstrated that due to the con- 
straint structure of the theory, the charge carrying holons 
need to have a hard-core interaction, which accounts for 
phase separation. The predicted compressibility and crit- 
ical doping are in accord with chemical potential shift 
measurements^'^'. These phase separation tendencies of 
the SU{2) gauge theory puts the door ajar for more in- 
tricate phenomena, like the stripe order-i^ii^ii^ that has 
been observed in experimentsiiii^iiSii^SiSi. 

In turn, this constraint structure turns out to be re- 
sponsible for an unexpected surprise: the superconduct- 
ing order parameter needs to be d -|- s-wave, instead of 
d. This interferes in an interesting way with the empiri- 
cal developments in high- Tc superconductivity. There is 
strong experimental evidence from c-axis tunneling22,f2^ 
and Raman scattering2ii2^ that in Bi2212 there is an 
s-wave component in the gap, which is in largeness 
comparable'^'^ to our prediction, and grows with doping^!, 
also in accord with our prediction. Further, 7r-phase shift 
experiments for YBCO point out that the s-wave compo- 
nent therein cannot be fully explained by the orthorom- 
bicity of the crystal^i^l. As far as we are aware, SU (2) 
gauge theory, in our formulation, is the only theory ex- 
plaining these results in at least an elegant way, and is 
demonstrated to be rooted in the way SU{2) gauge the- 
ories describe doping of holes. 

It appears that the above experimental findings are 
largely ignored because all existing mechanism theories 
predict either a d-wave or an s-wave, and the SU (2) 
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gauge theory is stand-only with regard to its insistence 
on a d + s-symmetry. In the narrow context of slave- 
Uke theories, Ogata and coworkers excluded d + s in 
the related context of Gutzwiller projected wave func- 
tion Ansatzes2Si2^. 

The organisation of this paper is as follows. In Section 
im we review the SU (2) gauge theory formulation of Wen 
and Lee, and mention the projective symmetry ideas of 
Wen. In Section lllll we take a somewhat warped perspec- 
tive on doping: to fully understand the effects of doping, 
we consider the empty limit, i.e., the doping x — 1. As 
academic as it might seem, this exercise is very instruc- 
tive as to the fact that the slave holons should have hard 
cores, and that the s-wave order parameter is induced by 
doping. It turns out that the s-wave component is equal 
to half times the density of dopant holons. We also show 
that the mean field theory is able to get both the energy 
and the single-electron propagator right. Furthermore, 
the mean field theory in the empty limit gives the inspi- 
ration for the mean field wave function to be exploited 
in our analysis for intermediate dopings in the following 
section. 

In that Section ITVl we derive a mean field free energy 
functional for non-zero doping in the grand canonical en- 
semble, to be able to account for phase separation. The 
mean field phase diagram is calculated. The phase sepa- 
ration properties of the hard-core holon SU (2) slave the- 
ory is further quantified by comparison of the compress- 
ibility with experiments on the chemical potential shift^. 
The agreement turns out to be very well. Also the node 
shift, caused by the s-wave admixture, is determined. 
The hard-core nature of the bosons leads to interesting 
and experimental falsifiable properties. Namely, the hole 
condensate leads to phonon modes in the incoherent part. 

In the last Section |V] the idea of the isospin spiral is 
introduced and quantified. We show that the mean field 
energy of the inhomogeneous state is only a bit above 
the energy of the homogeneous state. We argue that this 
should be viewed as an artefact of the t — J-model, and 
not of the isospin spiral state. We also show that the 
nodal fermions are not that much affected by the charge 
inhomogeneities as expected. In effect, the gap is very 
small, leading to only very small Umklapp scattering at 
the Fermi pockets. 

These interesting results form a promising motivation 
to study SU{2) gauge theory in more realistic models 
for the cuprates, the more so since it seems to be the 
only theory on the market predicting a, d + s-form of the 
superconducting order parameter. In fact, SU{2) gauge 
theory forms the bridge between, on the one hand, nodal 
fermions and spin liquid states, and on the other hand, 
striped superconductors. 



II. 5'C/(2)-SLAVE BOSON FORMULATION OF 
THE t J-MODEL 

As is widely accepted, the problem of high-temperature 
superconductivity, is the problem of doping a Mott insu- 
lator. The parent compounds of high- Tc are insulating, 
due to a large Coulomb repulsion. By removing electrons 
however , the charges get mobile, and it is believed that 
this physics is at the origin of the superconductivity. This 
forms the motivation to include the hopping of projected 
electrons Cia in the original Heisenberg Hamiltonian: 

Ht-j = Yl •^i^i' ^ l^^njJ-J2*'^ (cLcjT + /i.e.) 

<ij> ^ ^ ij 

(1) 

where the hoppings tij are the wave function overlaps 
of electrons at sites i and j. Without loss of generality, 
we take tij — t for nearest neighbour sites, and tij = 
otherwise. 

The idea of the SU{2) slave boson formulation, due 
to Wen and Lee^'^, is rooted in the observation by Af- 
fleck and Marston^^ that the spin system at half-filling 
in a fermionic 'spinon' representation is characterised by 
both the usual 'stay at home' U{1) gauge symmetry, and 
a local conjugation symmetry meaning that one can de- 
scribe the spin system equally well in terms of spinon 
particles and antiparticles. This idea is encoded in in- 
troducing the SU{2) doublet composed from the spinon 
operators /i| and 

In this language, the spin operators are given by 

S+ = li4^4^-4^4^) (3) 

s.^ = ^(V'lV'.-i). (4) 

Importantly, the Hilbert space of the t ~ J Hamilto- 
nian is formed from three states only: a spin up electron, 
c j I vac) , a spin down electron cj | vac) and a vacancy | vac) . 
Consequently, the electron operators in Eq.(IT]) are pro- 
jected electron operators — Cia{l ~ riia), where a 
denotes a spin opposite from a. The tilde is dropped 
from now on and the projection is kept implicit. Bear 
in mind that one should take care that all physics takes 
place in this projected Hilbert space! 

Let us now describe these electrons in the SU (2) gauge 
theory. Since in everyday life only the physical electrons 
are encountered, the electrons should be SU{2) singlets 
To construct these singlets, introduce the SU (2) doublet 
describing holons: 
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Then the appropriate SU{2) singlet describing the pro- 
jected electron is 

c^^ = -^hl^.^^iblj^. + byi), (6) 



V2 



V2' 



(7) 



The equations ^ as they stand, however, are not op- 
erator equalities, since the Hilbert space of the SU{2) 
theory is larger. To arrive at the correct physics, we 
have to impose constraint to make the mapping to the 
states of the original t — J-Hamiltonian exact, so that 
the Hilbert spaces are equal. This is achieved as follows. 
Since the electrons should be SU{2) singlets, we should 
require that the physical states |phys) of the SU{2) slave 
boson model obeyi 



Iphys) = 0. 



(8) 



There are precisely three states satisfying Eq. (|5]) . The 
first two states satisfying those, are fj- |vac) and |vac), 
corresponding the projected up and down electron in the 
t — J model. The description of the hole leads inevitably 
to the introduction of the holon doublet Eq.®. This is 
seen by the fact that the empty state |vac) is conversed 
into a spinon doublet by the projection operator in Eq. 
([5]), and vice versa. Henceforth, when an empty site is 
accompanied with a bi boson, and a doubly occupied site 
with a 62 boson, the hole in the t— J model is represented 
by 



|0), 



V2 



(9) 



It is easy to check that Eq. ([9]) does satisfy all the con- 
straints Eq. ((HI). Observe that we needed all the three 
constraints to arrive at this expression. In the original 
Wen and Lee formulation, the consequences of this fact 
have not been taken at face value. It will be demon- 
strated that this is not justified for higher dopings. 

In the expression for the hole, SU (2) gauge theory al- 
ready reveals some of its powers: it captures the fact as- 
sociated with the particle-hole symmetry intrinsic to spin 
that the singletness of pure vacuum should be treated on 
precisely the same footing as the spin singletness of either 
the empty or doubly occupied spinon configuration. 

In this way, it is shown how one can make a mapping 
from the slave boson operator states to the Hilbert space 
of the t — J model, by including exact constraints Eq.®. 
Solving for exact constraints is however extremely diffi- 
cult. 

In order to make progress, we are going to put forward 
a mean field theory, and treat the constraints on a mean 
field level. 

one may assume the existence of the fermionic vacuum 
expectation values 



X^j = if If]] + fljji) 



(10) 

(11) 



The first one is a hopping amplitude, inspired by the 
staggered flux spin liquid states proposed by Affleck and 
Marston^iii. The order parameter A is going to play the 
role of a superconducting amplitude with d-wave symme- 
try. 

Let us now assume that the expectation values Eq. 
(fTU)) and Eq. pT|) exist, and let us presuppose deconfine- 
ment, by neglecting ffuctuations in the gauge fields Aqj. 
This is equivalent to replacing of the exact constraints 
Eq. ([5]) by the mean field constraints 



(12) 



Let us first obtain a manifestly invariant SU{2) gauge 
invariant mean field theory, by grouping the mean fields 
as follows: 



A*, X^J 



(13) 



To actually calculate matters, we have to derive the 
slave boson version of the Hamiltonian Eq.([T]), with the 
decomposition Eq.®. To decouple terms quartic in the 
spinons, we use the spin liquid Ansatzes Eqns. pil)) and 
PT|) . In order to impose the three constraints Eq. ipg)) . 
we need to incorporate the Lagrange multipliers ag^ into 
the mean field Hamiltonian. Bearing these remarks in 
mind, it is a straightforward exercise to derive the mean 
field Hamiltonian 



mf 



3 J 



h.c. 



(14) 



where Uij was already defined in Eq. p^ . 

Of course, the Hamiltonian Eq. (fT4)) is manifestly 
SU(2) invariant under the transformation 



I I t 
g^a^QT g]. 



(15) 



Inspired by the spin liquid ideas of many people^i^i^, 
an idea having some experimental support^ , we intro- 
duce three mean field states, namely the staggered fiux 
phase^iii, the d-density wave state^i^i and the d-wave 
superconductor. We recall their descriptions here. 

The d-wave superconductor is in the projective sym- 
metry group represented by 



dSC 

Ui,i+y 



-XT 







1,2 



^0, 



< 5i > = < 62 0. 



(16) 
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The dispersion for the fermions is (at least for homo- 
geneous Ooi) readily calculated to be 



Eh 



Xfc = — ^ (cos + cosky)x, 
Ak = — ^(cos kx — cos ky)A. 



(17) 



One should notice that the spinous have gapless Dirac 
dispersions at the points {kx,ky) = {±^Tr,±^n), i.e., 
at those points we have nodal fermions consistent with 
experiments^^-'^. 

The following two phases also support these Dirac 
quasiparticles. The first one is the staggered flux phase, 

SFP 

Uu+y = -xr' + ihYA, 



1,2.3 r\ 

an = 0, 



< 5i > 



< 62 >= 0. 



(18) 



The third phase has Fermi pockets, with nodes radially 
shifted from (^tt, ^tt): 

dSC with pockets 

C/i,i+x = -xr^-ii-YA, 

Uu+y = -XT^ + ii-YA, 

al'^ = 0, ai^ 7^ 

< 61 > 7^ ,< 62 >= 0. (19) 

We remind the reader that the dSC state is referred to 
that way, since in the Hamiltonian Eq. (I14|) that partic- 
ular Uij couples /|i with fii. 

The SFP is a spin liquid, describing spinous hopping 
around the plaquettes of the square copper-oxide lattice. 
It breaks translation symmetry, since the hopping fluxes 

$.op = J Yl Arg({/ii) = (20) 

plaqucttc 

show a bipartite staggered pattern. By a Fourier trans- 
formation to momentum space, the dispersion is read- 
ily calculated. It turns out to be identical to the dSC- 
dispersion 



Ek = ^{Xk-4r + (Ak-al)^, 



(21) 



The three mean field phases above describe nodal 
fermions, but at different points in fc-space. 

For zero doping, however, all the three Ansatzes 
become the same, supporting Dirac quasiparticles at 
(^TT, ^tt), since the Lagrange multipliers and boson den- 
sities vanish. This is surprising, since the SFP breaks 
translation symmetry, whereas the dSC does not. In the 



framework of classical Landau-Ginzburg- Wilson theory, 
this is impossible: in general, different symmetry broken 
states give rise to different excitations. What is going 
on here? In fact, the two above states are two sides of 
the same coin. This is seen after applying the following 
site-dependent transformation. 



g^ = exp {-iji-Y'^^ 



(22) 



by which the SFP is mapped to the dSC. (The quantity / 
is defined as / = ix+iy) Put differently, the translational 
symmetry breaking of the SFP is just a gauge artefact. 

Due to the idea of WeniS, this is an expression of 
the fact that for zero doping, these states are members 
of the same projective symmetry group (PSG). This 
means that different mean field states describing the same 
physics, are connected to each other by gauge transfor- 
mations. Wen then declares them to have the same quan- 
tum order, a novel concept going beyond the standard 
Ginzburg-Landau- Wilson paradigm. 

In fact, all states connected to the dSC by a gauge 
transformation 



9^ 



exp I 



(23) 



are equivalent. This can be pictured nicely by the con- 
cept of the isospin sphere. An S'[/(2)-gauge group ele- 
ment gi can be written as follows: 



Zil -Zi2 
Zi2 <i 



(24) 



where the complex numbers zi are parametrised by three 
angles, viz., 

Zii = e*"'e *2 cos — , Zi2 = e'"'e' 2 sin—. (25) 

The z's are grouped in the vector Zi ~ {zn, 212). 

The isospin vector li turns out to be a useful definition: 

li zlrzi — (cos 0i sin 9i, sin sin 9i, cos 6i) (26) 

The angle can then be interpreted as the latitude on 
the isospin sphere, whereas the angle is the longitude, 
cf. Figure [III The north and south pole of the sphere 
correspond to a staggered flux phase, with A — B and 
B — A staggering respectively, while the equator corre- 
sponds to the c?-wave superconductor. For half filling, 
the rotations on the isospin sphere correspond to pure 
gauge transformations, meaning that spinon fiux phases 
and d-wave superconductors are gauge-equivalent. 

The way the SU{2) mean field theory is set up, is as 
follows. The spinous and holons are considered to be 
separate systems. As long as one rotates the spinous and 
holons together, SU{2) gauge symmetry gives the same 
mean field properties. If one fixes a gauge for the spinous, 
and then starts to rotate the holons independently, the 
mean field results for the energy will be different. The 
strategy we choose is to fix the spinon gauge at the dSC 
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A-B Staggered Flux 



Isospin vector 



d-wave SC 



simple, 



B-A Staggered Flux 



FIG. 1: The isospin sphere. States on the north and south 
pole are staggered flux phases, whereas states on the equator 
are dSC. In between, we have a DDW state. 



mean field state, whereas the holons will be rotated by 
the group element gi . It was shown in Eq. (j24p that gi can 
be decomposed into "Euler angles" ai^Qi and These 
can be encoded in the useful concept of isospin, as defined 
in (pS)) . We pictured this concept in the isospin sphere, cf. 
Fig. ini Only the angle Qi will turn out to be physical, 
whereas the other ones are gauge. For equator states, 
we have d = ^tt, corresponding to < bi >=< 62 >, 
whereas for 6^ = 0, < 61 >^ and < 62 >= 0. The 
latter means that the symmetry between empty sites and 
doubly occupied sites is broken. 

The concept of isospin will turn out to be important for 
the last section, since it can encode for inhomogeneous 
states as 6i might vary with lattice site i. But before 
turning to that topic, we need to cross some terrain. 



III. LESSONS FROM THE EMPTY LIMIT 

In this section, we investigate the consequences of the 
peculiar form of the expression for the hole in the SU (2) 
formulation. We will demonstrate in the first subsection 
that the spinon pair involved will give rise to an s wave 
component in the spinon pairing parameter. Further- 
more, we will show that the constraint equations lead 
to requirement of the bosons to have a hard core. This 
property will also take care of the mean field energy to 
be correct in the empty limit. In the last subsection it is 
demonstrated that the proposed mean field wave function 
even renders the single-electron propagator correctly. 

Let us first discuss the hard-core nature of the holons. 
The exact wave function describing the empty limit is 



|0)-n(^li+4/iV|.)|vac). (27) 



Deconfinement or spin-charge separation implies that the 
system loses its knowledge about the three particle cor- 
relation &i2/|i/ti ^^^^ ^'^^ '^^^ ^'^ ^ 
holon-spinon product wave function. The best choice is 
obviously 



vac) 



(28) 



This wave function still has to satisfy the mean-field 
version of the constraint Eq. ([5]) , 



(^T^. + ZiIr'/i,) =0, 



(29) 



where the brackets in this case stand for the expectation 
value relative to the state \0)mf- This brings us to the 
main point: the constraints are only satisfied when the 
bosons have hard cores. Indeed, the mean field wave 
function Eq. (|28p obeys 



(oi i^y^j, io) 



MF ■ 



h^T^h, |0) 



MF 



+ = 0. (30) 



If we were to take soft-core bosons, the state 
vac) would be possible. It does not sat- 



((^1)^ + 4) 



isfy the constraints, however, since then 

(0| Vlr^T^, |0)^,^ + (0| h,T^h, \0)j^,p =0+1^0. (31) 

So, if we are to take the Hilbert space constraints seri- 
ously, we need to accept that the holons have an infi- 
nite hard core. This is consistent with the fact that the 
holons carry electric charge. Since the Coulomb repulsion 
is taken infinite to arrive at the t — J model in the first 
place, this means that there can be at most one 6-boson 
per site. 

The reason to stress this hard-core nature of the 
bosons, is that in the original formulation^, the bosons 
were taken to be non-interacting. The empty limit ex- 
ercise, being transparent in this regard, shows that this 
is inconsistent for appreciable dopings. This is also il- 
lustrated by the fact that the hard-core bosons give the 
correct energy in the empty limit. Indeed, let us first 
show that with the original projected electrons, the ex- 
act energy of the t — J Hamiltonian Eq. ^ is zero in 
this limit. The empty state is described in the projected 
electron formulation simply by |vac), implying vanishing 
rii on this state. Since the spin operator S' is given by 
'^\a\'^ai3^ip^ the energy of the Heisenberg part vanishes. 
Since there are no electrons around, the hopping part 
vanishes as well, and we conclude that the total energy 
vanishes. 

Let us now demonstrate that the Ansatz Eq. (pS)) yields 
the same result. Firstly, the Heisenberg term vanishes. 
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The only component of spin that could contribute is the 
I — 3 component, since the others vanish when acting 

on ^1 + |vac) . However, since is the number 

of up spins minus the number of down spins, it vanishes 
as well. Furthermore, since the number operator ni in 
the slave boson representation reads c|^CiQ, = ^iplhihlipi, 
it also gives zero contribution, because of the hard core 
condition. Similarly, the hoppings vanish for the same 
reason. 

This would not be the case if the bosons were assumed 
to have no hard core, since for weakly interacting bosons 
there is no restriction on the hoppings. This is inconsis- 
tent with the t — J model, as hopping is only allowed be- 
tween occupied and empty sites. In conclusion, the hard 
core condition is a sufficient condition for the empty state 
to have the correct energy in the empty limit. 

So far, we have derived an exact expression for the hole 
creation operator in the SU{2) gauge theory, taking se- 
riously all three constraints. We considered the empty 
limit next, since it is easy to construct the mean field 
theory in this case. Treating the constraints correctly, 
we arrived at the conclusion that the bosons should be 
treated as having a hard core. In the next section, we will 
show how our mean field Ansatz Ea. (P5|) generalises to 
mean field wave functions describing intermediate dop- 
ings. 

A. Doping induces an s-wave order parameter 

In the previous section, we showed that a correct de- 
scription of the empty limit requires that the bosons 
should be treated as hard-core. The empty limit con- 
siderations leading us to that conclusion, turns out to 
be very useful to find out the structure of the mean field 
wave function at intermediate dopings. In fact, hard-core 
bosons are just like XY spins, and the straightforward 
generalisation of Eq. (|28|) becomes obvious, 

\'^o)hoions = n + l^^i'^Mi + vzbl^ij |vac) , (32) 

i 

where the complex numbers a and /3 obey the normali- 
sation condition jap -|- = 1. 

To already harvest some results from our considera- 
tions, let us consider the saddle point Lagrange multiplier 
equations 

/^(i/„/)=0, / = 1,2,3. (33) 

These are precisely the mean field constraint equations 
Eq.dMl): 

(./tViU/i./T^) = {b\M + blbu) (34) 
{fhfl ~ M^^) - {bhb2r ~ blbu) (35) 
(fLU-l) = (blM-blbu). (36) 



These constraint equations already convey an important 
message. The third equation tells us something about the 
deviation from half-filling, which was an important point 
already made by Lee, Wen and Nagaosa? . As soon as the 
average fermion occupation number deviates from unity, 
i.e., deviates from half filling, there is a difference between 
bi bosons and 62 bosons. In plain physics language: as 
soon as Fermi pockets form, the difference between empty 
sites and spin pair singlets becomes physical. 

The first two equations acquire a novel interpretation. 
For non-zero dopings, the holon expectation values are 
non-zero as well. However, looking at the left-hand side 
of the equations, one needs to conclude that a super- 
fluid order parameter appears with an s-wave structure. 
In other words, taking seriously all constraint equations, 
and having convinced oneself that doping must be de- 
scribed by a superposition of both empty and doubly 
occupied sites, one has to face an extra order parameter 
with a superconducting s-wave symmetry. Rephrased in 
physical language: within the framework of SU{2) the- 
ory, doping induces s-wave pairing. 

One could argue that there are some left-over degrees 
of freedom, so that one could gauge away the s-wave 
component. This is not the case, however. Let us ex- 
ploit the isospin representation, introduced in Section 
ITTl Using the isospin angles (p and 6, cf. Fig. |lll we 
parametrise the holon wave function (j32p by Ui = cos(^) 
and Vi = sin(^). Further, choose Pi /^ic"''' such 
that f3i is real. This parametrisation is instructive, since 
Oi = ^ makes the expectation values for bi with vacan- 
cies indistinguishable from 62 with a spinon pair, repro- 
ducing the particle- hole symmetric empty state Ea. (p8| . 
Moreover, this corroborates the point that the equator 
on the S'C/(2)-isospin sphere {i.e.,9i = ^ ) corresponds to 
the particle-hole symmetric d-wave superconductor. Cal- 
culating the expectation values explicitly, the equations 
Eq.dSHI) become 

(fh-fh + fufv^) = |/3o.|'sin(0Ocos(^,) (37) 

{.flfl-fu.h) - |/3o.psin(0,)sin(^,) (38) 

= Mcos{e,). (39) 

On the one hand, this illustrates once again that equa- 
tor states are particle- hole symmetric, cf. Eg. (15^ . Sec- 
ondly, and more importantly, we see that there is no way 
to gauge away the s-wave component. In other words, 
as soon as there is a superconducting order parameter 
(sin(0) 0), there is an s-wave component linearly in- 
creasing with doping x: 

A,^^xsm{e). (40) 

The only freedom is to choose its phase to be real by 
choosing (pi = 0, implying zero ctpj^". This is a first result 
of our empty-limit exercise, which at first sight looks triv- 
ial. Conversely, the first equation tells us that we cannot 
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neglect the Lagrange multiplier Qqi, which accounts for 
the s-wave admixture. This has been ignored in the orig- 
inal formulations of the mean field theoryl>^>a. This flaw 
leads to a disaster, as we will show in the next subsection. 



B. The empty limit in mean- field theory 

One could wonder if it is really wrong to leave out the 
first constraint. Probably one remains very closely to the 
"correct" mean field state when ignoring it? This is not 
the case: it leads to nonsensical results. The bright side 
is the ease by which the constraint ajj is incorporated. 

Let us fix the holon density (h\hij — 1, and choose the 
Hubbard-Stratonovich fields to be homogeneous, Ay = 
^iXij = X- Since the empty state makes no distinction 
between empty and doubly occupied sites, we have 9i ~ 
= ^TT. Inserting these assumptions into Eq. p^ . we 
obtain an energy density functional £",„./ for the empty 
limit. 

Let us first ignore a\i. Then the mean field equations 
are 



multiplier is of central importance, and the mean fields 
vanish, as they should. Substituting this solution in the 
Hamiltonian Eg . ([T4|) . we recover the correct energy for 
the empty limit. In other words, things go dramatically 



wrong if aj is ignored. In Section IIVI we will show that 
our mean field theory is performing well as an energy 
density functional at intermediate dopings. To calculate 
dynamic properties, one has to be careful in choosing the 
correct mean field approach. 



'mf 



dx 



dE, 



'mf 



dA 



= 2X = E 

k 

= 2A = E 



x(cos kx + cos ky)^ 
Eh 



A (cos kx — cos ky — aj) (cos k^ — cos 



Ek 



C. Dynamical properties of the empty limit 

It is interesting to consider the dynamic properties fol- 
lowing from the slave boson theory in the empty limit. 
It turns out that although this theory is a good energy 
density functional, it is less trustworthy with regard to 
dynamical properties revealed through the propagators. 
This is of course due to the mean field treatment ignor- 
ing fluctuations of the gauge fields. Indeed, the Lagrange 
multipliers should be given dynamics, causing fluctu- 
ations confinement of the spinous and holons to electrons. 
Since we ignored the fluctuations, we can not expect the 
mean field theory to describe the physical electron of the 
(^-'^^mpty limit. 

The starting point for the study of electron dynamics 
^l^ip the single electron propagator 



al = 0. 

where the dispersion iJj. is given by 



(43) 



G{x,y;t-t') 



{T{cx^(t)cl^{m (46) 
-.H(t-t')Act^)e(i_i') 



Eh 



^ ^ {Xk - alY + {Ak - al)\ 
Xk = — —{coskx + cos ky)x, 
Ak — — ^(cos kx — cos ky)A. 



-((c:.e 



Here the Cia- again describe projected electron operators. 
Since in the empty limit there are by definition no elec- 
trons in the vacuum, 



(44) 



G{x,y;t-t') = J2 



^ik' x — iky 



c/c'ie 



-iH(t-t')/hJ \ 

^k]l 



In the empty limit, the holons cannot move, so there are 
no mean field equations and dispersions governing those. 
The above mean field equations can be solved numerically 
to yield the unphysical result x = ^ — identical to 
the result for half filling. But this is clearly nonsense: the 
empty limit is neither a spin liquid nor a superconductor. 
Also, since A and x ^-re non-zero, the total energy will 
be nonzero, in flagrant contrast with the correct result 
being zero, as pointed out earlier. 

Taking aj into account, however, the above mean field 
equations are extended with the saddle point equation 
for Oq, 



kk' 



(47) 



Let us first show that for the exact expression ([77|) 
describing the empty limit of the t— J model, one obtains 
a free particle dispersion. 

We need to know the time evolution operator 



-iH{t~t')/h 



The Hamiltonian operator has no Heisen- 



berg part for projected electrons, and it also vanishes 
on the state c^^jlO), where |0) is the wave function 
Ea. (P7|) . So we only need to calculate the effect of 
Ht = Eij ^jhih^j'ipj on 



1 = 



E 



(A(cos kx — cos ky) 
E'k 



(45) 



I empty) 



^ \ irnx ( 



(&lm/L + 62,n/irn)|0). (48) 



where the number 1 is the boson density. Solving the new 
system of equations numerically, we obtain the correct 
result X = A = and a\ = 



In Fourier space, the result is 



Therefore, the Lagrange {ck'^Htcj,^) = —2tSkk'{coskx +cosky) = SkSkk' 



(49) 
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Including the chemical potential and using a contour in- 
tegral expression for Q{t — t'), we conclude that the prop- 
agator is 



hw - {sk - Ai) + 17] ' 



(50) 



the correct result for the propagator of a free particle. 
This means that the wave function Eq. ((27| encodes the 
right physics, as expected. 

What performance can be expected from the mean field 
theory when asked dynamical questions? Let us first con- 
sider the case for the mean field Haniiltonian Eq. 
Without a", one obtains the expected errors, namely 
a d-wave dispersion. Repairing this with a nonzero 
leads, however, to just partially good news. We showed 
that the d-wave dispersion vanishes, which is physical. 
The bad news, however, is that both x and A are zero. 
This absence of kinematics would lead to a dispersionless 
spinon spectrum. This is not what one would expect, 
since shooting electrons into the void should behave as a 
free electron cosine dispersion. 

This motivates the second approach: let us apply only 
mean-field theory at the wavefunction level, without in- 
troducing the Hubbard-Stratonovich fields Eq. lTUl) and 
Ea. ([TT|) . This means that the Hamiltonian then reads 



<rj> 



In the empty limit, we just have one contribution to 
the propagator 



Gmf{x,y;t-t') = ^e*^-'^- 



iky 



-iH{t~t')/hj 



kk' 



xQ{t-t') 



(52) 



The expectation values are taken with respect to the 
empty-limit mean field state Eq. ((28)) . denoted by \0)^,.[p. 
It is convenient to first calculate how the electron opera- 
tor acts on the empty limit mean field state, 

cIt|0)mf = -^{bixfl^ +bi.J^i) 



{h\,T^hx) zero. Furthermore, the fermionic expectation 
values vanish as well. As an example, 



\Plfy^fxiPxmi 



{^\PlfyJxiPx,y\Q)MF^^ 



(54) 



both for X = y and x ^ y. The other Vs are checked 
similarly. This implies that the constraint terms in the 
Hamiltonian Eq. (|5ip vanish. 

Then we demonstrate that the Heisenberg term also 
vanishes on c^^ \Q)j^.jp. This can be seen easily since Si ■ 

Sjfj-^Pj |vac) = Sjfj-^-SiPj |vac) = 0, by the same reason 
why the Heisenberg term of the mean field energy on the 
empty state vanishes. 

Only the hopping terms are non-trivial. This is intu- 
itively clear, since shooting an electron in the empty sam- 
ple, removes one holon, allowing the rest to move. This 
is corroborated by calculating the matrix elements of the 
exact hopping Hamiltonian Ht — ^l^ikjil^j: 

= (vac| Pl{b,y + b,yfy^fyMx + bl f^fl^) P X | vac) 
4 ^ 



;vac|(l + /,T/xi)(,44)|vac) 



/ , ^ixSjy, 

iv) 



(55) 



meaning that the electron can only hop to a nearest 
neighbour site if that site is empty. The Fourier trans- 
form then reads 



(0| Ck'THtcl^ \0)mf = -^tSk'k{coskx + cosky) = SkSk'k, 

(56) 

the free particle dispersion. This leads to the same free 
particle propagator as in the exact case. 



Gmf{k,L0) 



1 



huj — {sk — fJ') + if] 



(57) 



In other words, when one treats the slave boson Hamil- 
tonian exactly, the mean field wave function for the 
empty limit gives the correct single-electron propagator. 



IV. THE PHASE SEPARATED d + s-WAVE 
SUPERCONDUCTOR 



-^fx]P- |vac) 



(53) 



with the definition = I],^, Wu + + fljh)- 

First we demonstrate that the state c^^| |0)m_f ^ 
physical state, i.e., it satisfies the constraint Eg. p^ . It 
does, since baxPx |vac) = makes the expectation value 



In the previous section, we have discussed how one 
should describe doping within the context of 51/(2) gauge 
theory. By studying this carefully, we convinced our- 
selves of two important lessons. The first one is that a 
hole in the t — J model is described in the SU (2) theory 
by a superposition of a vacancy and a spin pair singlet 
state. This implied that doping induces s-wave pairing. 
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The second message is that the holons describing dop- 
ing should be hard core bosons, instead of gaseous, non- 
interacting particles employed by Wen and Lee. This 
hard-core is necessary condition in order to account for 
the fact that there can be at most one charge per site, as 
has been clear from very the beginning. 

Our results are summarised in the mean-field phase 
diagram, reflecting the phase separated d + s-wave su- 
perconductor. Another ramification we make quantita- 
tive. The s-wave admixture in the superconducting gap is 
shown to shift the gap nodes along the Fermi surface. We 
predict how this node shift behaves as a function of dop- 
ing, an effect which might be just within the resolution of 
present day angle resolved photo-emission experiments. 

At the end of this section, we will discuss the actual 
meaning of the mean field states of SU (2) mean field the- 
ory, with respect to the question which mean field states 
are superconductors, and which are not. The distinction 
will be made by the absence or presence of the Meissner 
effect. 



A. SU{2) energy density functional 

In the previous sections, we already set the stage for 
the mean field theory description of the doped Mott in- 
sulator. Now we are in the position to derive the en- 
ergy density functional. Inspired by the empty limit, our 
starting point is the mean field wave function, 

|*o)mf = n + + |vac) |F) . 

(58) 

The ket \F) describes the many body spinon state, and 
the boson density |/3p is the density of physical holes. 

The important point of SU (2) gauge theory is that the 
particle-hole symmetry is broken upon doping. Indeed, a 
hole is described by a superposition of vacancies and spin 
pair singlets accompanied by their own boson. In 

the particle-hole symmetric state, the bi and &2 boson 
are equal, meaning that this should correspond with 9 — 
^TT. This is the motivation for the parametrisation Ui = 
coa{^9i) and Vi — sin{^9i). 

The phase f is the same as the phase tpi of gi, which 
is gauge. From now on, we gauge cpi = everywhere. 
The transformation ([22|l mapping the SFP into the dSC 



corresponds with ipi — -^tt + In, as the reader can verify. 

For theoreticians, it is natural to first study spatially 
homogeneous mean field states. However, in the course 
of time it has become clear that strongly interacting elec- 
tron systems tend to form inhomogeneous states, like 
stripes. Being aware of this complication, let us neverthe- 
less study homogeneous states. Still, this exercise turns 
out to be instructive, in this regard. The reason is that 
we treat the Hamiltonian Eq. (|14p in the grand canon- 
ical ensemble, instead of the canonical ensemble. The 
original formulation of the SU{2) mean field theorji^i^ 
rested on the canonical ensemble as well. To account for 



condensation of the holons at finite doping, the temper- 
ature was taken to be finite, to find out that the particle 
number constraint leads to Bose- Einstein condensation of 
the holons, by treating /i simply as a Lagrange multiplier. 
We prefer a different approach, since considering first fi- 
nite temperature is a detour given in by the unphysical 
assumption htat the holons form a non-interacting gas. 
On the other hand, hard-core bosons are not only more 
physical, but also easy to treat at zero temperature. In 
Section IIIII we made the point that the bosons are in- 
teracting, making them superfiuid at zero temperature. 
The last motive is the possibility of phase separation, 
i.e., the possibility of coexistence of phases with different 
densities at the same chemical potential in the same vol- 
ume, giving rise to the need of performing the Maxwell 
construction. To anticipate this, it is necessary to take 
the chemical potential for the holons as control parame- 
ter, instead as the Lagrange multiplier enforcing a fixed 
density. 

Let us substitute the mean field wave function Eq. ([55]) 
in the slave Hamiltonian Ea. (jl4p . we obtain an expression 
for the mean field energy per site {Hmj) — ^mf- 
is the number of lattice sites.) 



N ^ " 



k 



lAp) 



(59) 



- 2tx\a\^W - it^ + ao smO + al cos9)\P\^. 

The homogeneity of the Ansatz is expressed in the fact 
that the lattice site subscript i is dropped. It is impor- 
tant to observe that the isospin latitude angle appears in 
the mean-field energy, expressing the fact that for non- 
zero doping, SU{2) gauge symmetry is broken. Hence 9 
is not gauge, but has acquired physical meaning. The 
kinetic part of the energy is gauge invariant, leading to 
the fact that only the hopping amplitude shows up in the 
holon hoppings. From this density functional, we derive 
the saddle point equations for the dynamical variables 
X, A, ao,ag and the hole density = p{x) — ^ — 
To simplify matters a bit, we take the isospin angle 9 
as an external parameter, controlling the density of 62 
relative to bi. 

The saddle point equations in the grand canonical en- 
semble are easily derived, with the homogeneous forms 
of Eq.dSll) and Eq.(l39l), 

X.(x.-aB)^2fiiA,(^)(l-p(x)) 



2X 
2A 

P{x) sin(6') 
Pix) cos(6') 



Ek 

Afc(Afc -4) 
Ek 



3Jdx 



E 



(Afc - ah) 
Ek 



iXk - flo) 



= 



Eh 



, , , , , 1 / /i -I- ai sin 6* -I- On cos ( 
Pix) ( Pix) - 2 (^1 + ^ 
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As already discussed after Eq. p7|) . these equations 
give rise to an important law which provides a linear 
relationship between the s-wave spinon pairing As = 

■filfh) ttis doping X ^ pix), 

A, = ixsin(0), (60) 

as a direct ramification of the full constraint structure. 

In order to find the solutions to the saddle point equa- 
tions Eq. ()60p . the energy Eq. (|60|) is minimised numeri- 
cally using the simulated annealing mcthod'^^-'^'* . 

We point out that the fifth equation admits both zero 
and non-zero solutions for the density p. As a function 
of /Lt, the mean field energy will tell which one is more 
favourable. In Section llVBl we will show that the system 
chooses between these two by a first order phase transi- 
tion, and not a second order one! This means that the 
mean field theory [60] implies a phase separation regime, 
for the usual Maxwell construction reasons, when trans- 
forming to the canonical ensemble. 

B. The SU{2) mean field phase diagram 

We have now arrived at the point where we can col- 
lect the results. Our first step was to prove that one 
needs all SU{2) constraints to project onto the t — J 
model Hilbert space, while this constraint structure is 
also required for the mean field description of the SU (2) 
gauge theory. This will bring us to the first result: that 
the superconducting order parameter needs to have an 
s-wave component! Then we spent effort in proving that 
the holons need to be hard-core in order to respect the 
full set of SU{2) constraints. We now show that the su- 
perfluid hard-core holon condensate displays phase sep- 
aration behaviour, as expected for hard-core interacting 
systems. We thereby achieve an intrinsic connection be- 
tween slave theories on the one hand, and the observation 
of inhomogeneous states on the other hand, a connection 
that is traditionally considered as absent. 

Let us first discuss the numerical verification of a claim 
inferred in the literature^, namely that at finite hole den- 
sity the superconducting state {9 — tt/2) is preferred 
over the flux phase {9 = 0, tt) (cf. the inset in Fig. 1). 
This is a natural ramification of the breaking of SU (2)- 
symmetry for non-zero doping. The 9 ^ ^ states, char- 
acterised by Oq = 0, are energetically more favourable, 
being consistent with the instability of the flux state to- 
wards d-wave superconductivity, as already understood 
in the early nineties^^. 

The reader might already have noticed a kink in cmf 
as a function of doping in the inset figure HVBl In other 
words: there is a first-order phase transition. Indeed, our 
main result is that generically this mean-field theory pre- 
dicts phase separation at small chemical potential. The 
system stays initially at half-filling and pending the ratio 
of J/t at some finite p, a level crossing takes place to a 
state with a finite doping level, cf. Fig.l. 




■10 -0 8 -0.6 
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doping 



FIG. 2: Chemical potential shift as a function of doping, 
showing the phase separation behaviour of the reformulated 
mean field theory. Indeed, the chemical potential starts to 
shift for appreciable dopings only. The blue line (colour on- 
line) are the numerical results for J/t = 0.1, and agrees very 
well with the experimental results from Fujimorii^ (dotted 
line). The red line are the results for J/t — 0.3. The critical 
doping changes, but not the compressibility. The inset shows 
that ignoring aj, i.e., ignoring the s-wave component, gives 
a false vacuum. Indeed, there is a positive relative energy 
difference between the pure dSC and our mean field theory, 
growing with doping (lower line). 



We stress here that this phase separation behaviour 
is eventually coming from the hard-core nature of the 
holons: also the "wrong" mean field states (SFP and pure 
dSC with pockets) exhibit first order behaviour. We con- 
clude that, due to the Hilbert space restrictions on the 
SU{2) description of the holons, the theory insists on 
inhomogeneous states for low doping. As a function of 
increasing J/t the width of this phase separation regime 
is increasing (see Fig. 2) and we find that for J/t ~ 4 
the phase separation is complete. This is consistent with 
exact diagonalisation studies on the t-J model indicat- 
ing a complete phase separation for J/t > 3.5^^^. This 
is quite remarkable and it reveals that the gauge mean 
field theory has to be a remarkably accurate quantitative 
theory of the density functional kind: it is a good de- 
scription of the empty limit and the Mott insulator, and 
gives a fair prediction of phase separation tendencies. We 
stress that as a rule less severe demands on physical re- 
ality are put on density functional theory , instead of full 
dynamical theories. 

The significance of this finding is that for this most 
sophisticated version of spin-charge separation theories, 
phase separation is natural feature, as it is in the em- 
pirical reality. It is well understood that these macro- 
scopic phase separated states are an artefact of the 
oversimplified t — J model. By taking the long-range 
Coulomb interaction into account this will turn immedi- 
ately into the microscopic inhomogeneityi^r^f^, of the 
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2.0 3.0 4.0 J/t 



FIG. 3: Phase diagram as a function of doping x and the 
ratio J/t, according to the Maxwell construction. For dopings 
below a critical doping line Xc{J/t), homogeneous states are 
metastable against phase separation, stripes etc. For J ~ At, 
total phase separation takes place and the system becomes a 
mixture of Mott insulating and empty regions. 




kind that are seen in STM-experiments-^ . To see how 
well this slave theory handles the 'big numbers' in this 
regard, we show in Fig. 1 the electronic incompressibil- 
ity 1/k = = If according to the SU{2) theory, 

to find that it compares remarkably well with the exper- 
imental results due to Fujimori and coworkers^. Firstly, 
it is seen that independent of the ratio J/t, the compress- 
ibility is right on spot of the experiments: the slope of /i 
vs. holon density is the same as the measured slope. The 
doping at which phase separation occurs, namely 13 %, 
is correct only for the value J/t = 0.1, which is too small. 
Indeed, from ARPES measurements a ratio of J/t = 0.3 
is more realistic, but for those samples phase separation 
takes place at dopings of about 17 %, and not 21 %, as 
found in the SU{2) mean field theory for J/t = 0.3. 

Let us now focus on the nature of the superconduct- 
ing order parameter found elevated doping levels. As 
expected, the s-wave component becomes increasingly 
important, cf. Eq. (|60p . To further emphasize this , 
we compare in the inset of Fig. 1 the energy of a state 
where we have fixed the Lagrange multiplier aj = such 
that the s-wave component vanishes, with the best d + s 
mean-field state, finding that the former is indeed a false 
vacuum. To mimick the average behavior of the super- 
conducting order parameter also in the micro-phase sep- 
arated states at low dopings, we calculate matters now 
in the false (uniform) vacuum of the canonical ensem- 
ble, fixing the average density, simplifying the mean field 
equations Ea. (|60p . Indeed, p{x) becomes now a fixed p. 
Furthermore, since the state 9 = ^tt is favoured, we take 
the mean field Ansatz Uij to be the SC one, as stated be- 
fore, and consequently we have aj] = 0. This enables us 
to map out the phase diagram as a function of doping and 
the ratio J/t, leading to three phases. The first, for low 
doping, is the phase-separated, underdoped c? -I- s-wave su- 
perconductor. It is a mixture of charged, superconduct- 
ing islands in an insulating sea without charges, where 
the full SU{2) symmetry is restored. This reminds the 



FIG. 4: Zero temperature phase diagram of the t — J model 
according to 5'il7(2)-mean field theory. It incorporates three 
phases, viz. the phase separated d + s-wave superconductor 
for low dopings, the homogeneous d -\- s-wave superconductor 
for intermediate dopings, and a homogeneous s-wave super- 
conductor at high dopings. The bold line shows indicates the 
border of the phase separation region. The phase separation 
tendency grows for increasing J/t, to become complete at 
J/t = 4. The colors indicate the total superconducting gap. 
For zero doping, there is only a d-wave component, whereas 
the s-wave admixture grows linearly with doping, so that the 
total gap is non-zero even beyond the critical x — J/t line 
(rightmost line), where Ad vanishes. The inset shows the d- 
wave and the growth of the s-wave component separately for 
J/t = 0.3. The blue line indicates the doping level where 
phase separation terminates, computed by imposing unifor- 
mity (canonical ensemble) for x < 0.2. 

reader of the STM-pictures from S.C. Davis' group^. For 
intermediate dopings, the homogeneous, overdoped d+s- 
wave superconductor is found, whereas for high dopings, 
the d-wave gap vanishes, leaving behind a pure s-wave 
superconductor. Hence, although d-wave superconduc- 
tivity leads to an s-wave admixture, the reverse is not 
true. 

We find that the regime where phase separation is im- 
portant, the s-wave component is not negligible. This is 
consistent with Raman measurementa^i^ where the su- 
perconducting gap was found to have both d- and s-wave 
components. Although screening effects in Raman scat- 
tering make it difficult to compare our results directly to 
theirs, their results indicate that the ratio r = As/A^ 
grows with doping, as it does in our approach. Looking 
to Figure 3., we find in the phase separated overdoped 
regime s-wave admixtures of about r = 10 — 20%, con- 
sistent with c-axis tunneling experiments^^. 

We predict that at a doping level that appears to be 
higher than can be achieved in cuprate crystals a phase 
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FIG. 5: In this figure, we show how the s-wave admixture 
r = As /Ad grows as function of doping (red hue) for J — 
0.3t. The insets show the symmetry of the gap function in 
reciprocal space, indicating the shift of the gap node along the 
Fermi surface, coinciding with the Brillouin zone boundary Y- 
X. This shift is denoted by S, and grows with doping. The 
blue line plots 5 as function of doping in units of ir/a, where 
a is the lattice spacing. Up to a doping of 20% the results 
of the false vacuum homogeneous solution in the canonical 
ensemble are used, exploiting the Maxwell construction. To 
remind the reader of the fact that the homogeneous states are 
false vacuum states, the lines are dashed. Finally, the lines 
end where the d-wave gap vanishes. Here r is infinite and the 
Fermi surface is fully gapped. 

transition occurs to a pure s-wave superconductor. As 
we already alluded to, the gauge fluctuations should be- 
come more severe as well, for increasing doping and at 
some doping level a transition should occur to a confining 
"electron like" system. 

As we learned from the empty limit, it still make more 
sense than the result obtained by disregarding the first 
constraint equation Eq.(l37|). since that would lead to the 
unphysical result % = A = , giving the wrong vacuum 
energy, as we discussed in Section HIIl In other words, our 
mean field theory is a remarkably good density functional 
theory, but inevitably the theory fails completely in dy- 
namic regards. In the confined phase, at sufficiently high 
dopings, we need an approach which is qualitatively dif- 
ferent from slave theories, let alone that one can get away 
with the mean field version. On the other hand, in the 
superconducting doping regime, we find some promising 
experimental support for our results with regard to the 
d-\-s structure of the order parameter, meaning that con- 
finement physics might not be overwhelmingly important 
in the low-doping part of the phase diagram. 

Having said this, the experimental support for an s- 
wave admixture makes it possible to come up with a fal- 



sifiable prediction for photo-emission experiments. The 
s-wave component induces a shift in the nodes, as can 
be inferred from the spinon dispersion Ek, Eq. (fT7|) : the 
hopping vanishes along the line ky — t: — k^, so that 
the locus of the node can be readily calculated to be 
ky — ^TT ~ arccos(aQ/(3 J/2)A), which shows a doping- 
dependent behaviour as well, cf. the blue line in Figure 
3. The node shifts might be able to explain the U-shaped 
gap as measured in Bi2212^, since the twinning of sam- 
ples mixes regions of node shifts +S with —S, smearing 
out the V-shape of the gap to a U-shape. 

C. Single-electron dynamics 

It is also possible to interrogate the single-electron dy- 
namics of the superconducting state. Wen, Lee, Nagaosa 
and Ng have addressed this issue for non-interacting 
holons^. Since we showed that one cannot ignore the 
hard-core nature of the holons for appreciable dopings, 
it is interesting to investigate what the electron corre- 
lation function looks like. We will demonstrate that it 
consists of two parts, ubiquitous in hard-core interacting 
Bose systems-^: one is a coherent condensate piece, and 
the second is an incoherent piece displaying sound modes 
for long wavelengths. This novel result is entirely due to 
the hard-core nature of the bosons, and might lead to 
falsifiable experimental predictions. 

The calculation can be a bit simplified by the fact that 
in the superconducting state 9 = ^tt, hence we write 
6i = 62 = b. Then the electron propagator 

G{k,t~t')^{{ckit)cl{t'))) (61) 

is the convolution of bosonic operators bp and fk+p, which 
can be seen by Fourier transforming the expression for 
the electron — bl^fi^ + bl^fj^. The time resolved 
propagator G(fc, lo) can be calculated by using the spinon 
coherence factors Uk+p and Vk+p from the Bogoliubov di- 
agonalisation of the spinons 

/,T = UqJqO+Vq-fl^ (62) 
fll = -VqlqO + Wg7gl (63) 

Exploiting the Fourier transformation of the textbook 
result^^ for the boson propagator, the frequency resolved 
propagator then reads for finite repulsion U 



+ 



- \Uk+p - Vk+p\ 

{huj - hil + e{_|_p) + irj 

\Uk+p + ffc+pl 

{hw ~ hn + el^j^) - iri 



xoSip)S{n) 



h{Q + Qp 



h{n + ftp) 

^ ft2(f7p + f7)2 - {2Uxoy 



(64) 



Here, the fermion and boson dispersions are denoted by ej, and Hflp, respectively. Although the expression for 
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the propagator Eg. ([64)) seems lengthy at first, it con- 
veys a couple of important messages. The most obvi- 
ous one is the fact that the propagator is a convolution 
of bosons and fermions. The second one is the separa- 
tion between the condensate piece, proportional to xq, 
and an incoherent piece. The condensate piece leads to 
poles at the fermion dispersions, huj — sj^. The incoher- 
ent piece is even more interesting. Since for the bosons 
hflp = 2tx{coskx + cosfcj,), we have for the long wave- 
length limit and for fj, 2t (the requirement for a con- 
densate to exist) that the holon poles are located at 

hil ~ 4:tx\/2U^\k\ finite U. (65) 

The holon velocity of sound can be written in a form for 
hard-core bosons as well, since for finite U, x — ^{2tx + 
/i), whereas for hard core bosons 2tx + A* = 2x0. Hence 
the result Eq. (|65p can be generalised to infinite repulsion 
U, 

hn = 2{2tf^\k\ infinite U. (66) 

This means that for the incoherent piece of the propaga- 
tor Eq. ([64]) . the poles are located at 

huj = ei^p~v\k\, fc^O. (67) 

with the velocity of sound v = 2{2t)'^ ^/xq. 

In summary, we have argued that in order to achieve 
consistency in the SU (2) slave boson theory, one has to 
implement hard-core bosons instead of gaseous Bogoli- 
ubov bosons to describe the charge sector. As a result, 
one obtains phase separation at lower dopings consis- 
tent with the experimental observations. Also the com- 
pressibility matches very well. Furthermore, by inspect- 
ing the empty limit, we showed that an s-wave compo- 
nent in the superconducting order parameter is implied 
when d-wave superconductivity occurs, at least for dop- 
ings where homogeneous states exist, a ramification of 
the constraint structure. This finds its origin eventu- 
ally in the particle-hole symmetry central to the gauge 
structure of the SU (2) theory: to describe physical spin 
singlets, "no fermions" are indistinguishable from an "s- 
wave spinon pair" . This is reflected in the constraint 
equations, required to reduce the SU{2) Hilbert space 
to the Hilbert space of the t — Jmodel. The constraint 
equations tell us that as soon as d-wave superconductiv- 
ity emerges, one necessarily has an s-wave component. 
This s-wave admixture is in accord with Raman^i and c- 
axis tunneling experiments^!^. We also predict a node 
shift in the gap function, that might be measurable by 
photocmission. 



V. ISOSPIN SPIRAL STATES IN CUPRATES 

Having established the phase separation tendencies, 
new perspectives are opened as to what the role of stripes 



in the superconducting cuprates is. The understand- 
ing is that phase separation is a necessary condition for 
stripes to exist. On the other hand, phase separation 
has never been seen in the high-Tc's. In fact, there are 
other properties at work. The Mott insulator is an an- 
tiferromagnet, which makes it advantageous for holes to 
order in stripes instead of phase separated islands. In 
fact, stripes should be regarded as antiphase boundaries 
in the antiferromagnet^. 

In the past decades, many approaches to understand 
the Hubbard model in some slave boson representation 
are made. One is the large-^* expansion^i^, taking the 
limit of the spin value 5* ^ oo. The other is introduc- 
ing more than two flavours of spin, such that an SU{N) 
model is obtained . The large- iV limit leads to dimerised 
states^Si^, whereas the vacuum of the large-S" limit is 
the antiferromagnet^^^. The question arises if the group 
SU{2) is able to describe the "anti-phase boundariness" 
of the antiferromagnet, the more so since the antifer- 
romagnet and the dimerised state are incompressible, 
whereas the flux phases and the dSC spin liquids of the 
SU{2) theory are compressible. 

We propose a way in which the SU{2) mean field the- 
ory can describe anti-phase boundaries, within the spin 
liquid states descending from the Mott insulator. The 
first observation is that for zero doping, the dSC or SEP 
state is just a gauge fix within the same projective sym- 
metry group. Let us now consider a gauge in which the 
spinon mean field Uij Eq. (fT5|) rotates over the whole 
isospin sphere, 

U,, = exp{z{-iye,r') (^'^ exp(-z(-l)^a,ri) 

(68) 

by a harmonically varying isospin angle 

9i — Q ■ i — qix, ordering vector Q in a; direction. 

(69) 

In this way, the SEP state is smoothly connected to a 
dSC state. Observe that this state is in the same PSG 
for zero doping. Then the idea is that for underdoped 
samples this spiral state might be lower in energy than 
the phase separated state for the homogeneous d-l- s-wave 
superconductor. A cartoon representation is given in the 
figure [HI 

The peculiar feature of the isospin spiral is that in the 
SU{2) gauge theory the charge-density wave is made out 
of a superconductor, which is not the case in the large-S' 
limit. The antiferromagnetic domains are now replaced 
by a spin liquid, carrying nodal fermions, a feature ab- 
sent in the antiferromagnet. As the nodal fermions exist 
in both the dSC and SE phases, a very promising per- 
spective is opened up, supported by experiments. 

The first support comes from the results from Euji- 
mori and coworkers^ and Z.X. Shen and collaborators^ 
for LaSCO. The chemical potential shift measurements 
of Eujimori in underdoped LaSCO are compatible with 
the existence of charge order, with Shen finding sim- 
ilar results. Eurthermore, in the Nd-doped cuprates. 
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FIG. 6: A pictorial representation of the unit cell of the 
isospin spiral. The green areas correspond with superconduct- 
ing stripe regions (isospin angle Qi =0), blue is the AB flux 
phase {Q = and red the BA flux phase {6 — — ^tt). The 
drawn cos(2^i^) profile shows the boson density as a function 
of the x-coordinate in the unit cell. Here, di is the isospin 
angle of the fermionic mean field state Uij. Note that the 
bosonic isospin angle as defined in Section IIVI Eq. 1581 is in 
this case equal to Qhos = 9i — §7r. Hence, = (62) still holds 
in the superconducting state. It is seen that the hole-rich re- 
gion forms an antiphase boundary for the SF-liquid state in 
between the superconducting stripes. 



clear features of static stripes are measured already in 
the ninetiesi^ii^. On the other hand, existence of nodal 
fermions in underdoped cuprates is found as well^. The 
combination of these resuhs seem to indicate the coexis- 
tence of striped charge order with nodal fermions. This 
idea is backed by recent results from the group of J.C. 
Davisi, reporting that charge order and nodal fermions 
can coexist. 

The SU{2) gauge theory is able to capture 'Mottness', 
d-wave superconductivity and the protection of nodal 
fermions. The framework of the isospin spiral state in 
SU{2) mean field theory forms an excellent explanation 
to explain the mystery why nodal fermions should exist 
in a strongly correlated background. This is a promis- 
ing motivation to study the stability of the isospin spiral 
mean field states in underdoped cuprates. 

As a first attempt, we adapted our program such that 
it can incorporate inhomogeneous constraints, within the 
harmonic approximation. We used a 2x8 unit cell, as 
suggested in Figure Via simulated annealing, a self- 
consistent solution can be found for chemical potentials fi 
admitting inhomogeneous boson densities. In Figure ([7]), 
the mean field energy as a function of average doping is 
plotted , together with the energy of the homogeneous 
solution. 

It can be seen that the homogeneous states are bet- 
ter in energy in the homogeneous regime (x > 20% for 
J/t = 0.3). One expects that in the phase separation 
regime, the isospin spiral state wins. Unfortunately, the 
spiral state is energetically a little bit less favourable the 
homogeneous state, but this is an artefact of the unreal- 
istic t — J-model. Any model which is closer to the ex- 
perimental reality, admitting stripes in the underdoped 
regime, are expected to favour the spiral states. This 



-0.4 



FIG. 7: The mean-field energy of the isospin spiral state as 



a function of doping. The Q 



spiral state is taken, as 



depicted in Figure ([6}, such that the unit cell is 2 x 8. For the 
t — J-model with J — 0.3t, t — 0.44 meV, the homogeneous 
state is lower in energy, but not very much. Furthermore, one 
can see that the isospin spiral state does not allow nonzero 
doping in the phase separation regime. This should be viewed 
as an artefact of the t — J model, however: any more realistic 
model, admitting stripes, might make the isospin spiral state 
more favourable in energy in the phase separation regime. 



hope is corroborated by the fact that the energy differ- 
ences with the homogeneous state are small, i.e., a cou- 
ple of percent. The direction of future research is clear: 
do models supporting stripes, also support isospin spiral 
states in the SU (2) slave boson formulation? This ques- 
tion is of high importance in the underdoped regime. In 
that regime, it is interesting to see whether there can 
be made a "Yamada-plot" of the isospin speed Q as a 
hmction of doping x, since Q induces period 2(5-stripes. 

In addition, we can shed some light on the issue of 
the coexistence of nodal fermions with stripes. We deter- 
mined numerically the dispersion relations for the spinous 
for the Q — ^ spiral, on a 8 x 2 unit cell. We only took 
the spinon part of the propagator Eq. ([M)) into account. 
The result is shown in Fig.®. 

The surprising result is that there is no (appreciable) 
Umklapp scattering in the spinon sector, at this resolu- 
tion at least. It is interesting to calculate the full electron 
spectral density, taking both the inhomogeneities in the 
boson condensate and the incoherent phonon part into 
account. We expect that the Umklapp might vanish al- 
together. 

The Figure [5] might not be decisive, but very promising 
as to the idea that nodal fermions and stripes are not 
as mutually exclusive as thought for a long time in the 
high-Tc community. Further effort should be put in more 
precisely mapping out the nodal dispersion, especially in 
the underdoped regime. 

In summary, we discussed how the constraint struc- 
ture of the mean field SU{2) gauge theory of the doped 
Mott insulator forces the holons to have a hard core. 
This was demonstrated by considering the empty limit. 
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FIG. 8: The dispersion of the spinons in the Q ~ ^ spiral 
state, on a 8 X 2 unit cell. The dispersion is along the nodal 
direction. The energy resolution is Ahui/J = 1/160, and 
the momentum resolution is Ak/a = 1/128. The chemical 
potential was chosen such that the average doping was 26%, 
with variations (stripes) in hole density of 2%. The dispersion 
crosses the spinon Fermi surface at (|vr, iyr). It is seen that 
in contrast to usual belief, there is no appreciable Umklapp 
scattering. This result supports the idea that by the isospin 
spiral state, the nodal fermions coexist with stripes. 



The empty limit was shown to produce the correct mean 
field energy, using the hard-core condition, and the cor- 
rect single-electron propagator. As a surprise, we showed 
that an s-wave pairing in the spinon sector is inevitable, 
and grows linearly with doping. Then the mean field be- 
haviour for intermediate dopings was shown to capture 
phase separation tendencies in the underdoped regime, 
being quantitatively in accord with both numerical re- 
sults and experiments. The single-particle propagator 
was shown to display phonon modes. 

We put forward indications that the SU (2) gauge the- 
ory is able to capture both 'Mottness' and d-wave su- 
perconductivity in a unifying picture, introducing the 
concept of the isospin spiral. This idea might offer the 
possibility of the protection of nodal fermions, as the cal- 
culated spinon sector of the single electron spectral func- 
tion seems to indicate. The framework of the isospin spi- 
ral state in SU{2) mean field theory forms an excellent 
explanation to explain the mystery why nodal fermions 
should exist in a strongly correlated background. This 
is a promising motivation to study the stability of the 
isospin spiral mean field states in underdoped cuprates 
further, together with making predictions for ARPES 
measurements. 
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